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SUMMARY 

This is the third annual report on: Sea Ice-Atmosphere Interaction - Application of Multispectral 

Satellite Data in Polar Surface Energy Flux Estimates. 

The main emphasis during the past year was on: radiative flux estimates from satellite data; 

intercomparison of satellite and ground-based cloud amounts; radiative cloud forcing; calibration 

of AVHRR visible chann els and comparison of two satellite derived albedo data sets; and on flux 

modeling for leads. 

The highlights can be summarized as following: 

. Satellite cloud amounts are generally 5 - 35 % less than observed values over the entire Arctic. 

• Shortwave and longwave fluxes at the surface were modeled from satellite data on a monthly 
basis for the entire Arctic. Net longwave flux was found be negative throughout the year with 
a minimum value of -70 Wm'^ in the month of May. 

• Regional variations in radiative fluxes were discussed and related to variations in cloud 
amount. 

• Seasonal variations of shortwave-, longwave- and net-radiation cloud forcing at the surface 
were computed; positive net forcing values of 65 Wnr 2 were found for October through 
January, and the net cloud forcing was negative for the months May through August. 

• Using in situ radiosonde measurements narrow-band surface albedo values from AVHRR 
channels 1 and 2 can be estimated with an accuracy of 2 %, and 14 % respectively . A large 
discrepancy between modeled and observed radiance in AVHRR channel 2 is believed to result 
from a satellite sensor drift. 

• Albedo values as derived from two different satellite data sets show good agreement in spatial 
patterns for the entire Arctic. Albedo values derived from ISCCP were lower by 0.1 (10 %) 
compared to a manual DMSP visible band analysis. 

• Relationship between satellite foot prints and averaged turbulent fluxes was established for 
leads. The turbulent flux as derived from a 80 m field of view (FOV) decreases by 45 % for a 
satellite footprint of 640 m FOV. 

• Narrow leads have been found to be more effective in turbulent flux transfer than wide leads of 
the same aerial extent 


Officially the three year project has ended in April 1993. We have extended this project at no-cost 
through December 1993 to complete the analysis and parameterization of energy flux 
measurements over different ice thickness types. Sensible and latent heal fluxes were measured by 
K Steffen in Lancaster Sound during the two winters 1979 and 1981. This data set will now be 
analyzed by T. DeMaria (supervisor: K. Steffen). The results from this study will be crucial for 
the discussion of satellite-based turbulent flux estimates. So far no in situ flux measurement has 
been published for Arctic conditions throughout an entire winter (November though February). 
The final results of this study will be summarized and presented as a final report in early 1994. 
The main emphasis will be a comparison of different flux components for different ice thicknesses 
as compared to modeled results based on mean climatology. 
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1. ARCTIC CLOUDS AND RADIATION 

Jeffrey Key and Axel Schweiger 

1.1 ISCCP and Surface-Based Cloud Amounts 

The International Cloud Climatology Project (ISCCP) is now producing a global cloud product 
(Rossow et al., 1988). How valid are the reported quantities in the Arctic? This question is 
critical if we are to use this data product for studies of Arctic climate, for example in the 
estimation of radiative fluxes as described in the next section. As a first step in determining the 
validity of these data we have compared the satellite-based cloud amounts to surface-based 
observations. However, there are many questions concerning the validity of the surface-based 
cloud amounts as well, so it is not possible to draw conclusions in an absolute sense about either 
data set. Only a summary of the results is presented here; for more details see Schweiger and 
Key (in press). 

The comparison is made with monthly cloud statistics from an atlas of global cloud cover 
compiled by Warren et al. (1988; hereafter WAR). Cloud observations from ships for 1951-1981 
form the basis for the ocean area atlas. Since there is no direct overlap between the ISCCP and 
the WAR data sets, only monthly statistics were compared. Interannual variability of monthly 
cloudiness is small, in the range of 5% for July and up to 8% in January of the four ISCCP years 
1983 to 1990, so that errors due to non-overlapping climatologies should be of a lower 
magnitude than the other potential error sources. 

Figure 1 shows the annual variation of total cloud amounts in the ISCCP-C2 and WAR data sets 
for the region north of 62.5° latitude. Two separate cloud amounts are given for the ISCCP-C2 
set: the total cloud amount reported in the ISCCP-C2 data set and marginal cloud amount. The 
total cloud includes the marginal cloud. Marginal doud amounts are derived in the ISCCP 
procedure through the application of a different threshold and represent the first derivative of 
doud amount with respect to the threshold. These marginally cloudy pixds display radiance's 
that are near the dear sky radiance's and occur primarily in areas of broken doudiness and where 
the contrast between doudy and clear conditions is low. Such conditions prevail in the polar 
regions and the marginal cloud amount is presented as an indicator of the levd of uncertainty in 
the thresholds. As can be seen from Figure 1, this level of uncertainty is greater during winter. 

The satellite cloud amounts are generally 5 - 35% less than the surface observations over the 
entire Arctic. However, regional differences may be as high as 45%. During July the 
surface-based cloud amounts for the central Arctic are 40% greater than the satellite-based 
values. In winter the ISCCP climatology in the central Arctic agree to within 10% of the surface 
observations. While satellite-derived doudiness during summer seems to reflect the conservative 
approach of the ISCCP algorithm and probably represents an underestimate, there is considerable 
uncertainty with respect to the winter doudiness. Ice crystal precipitation could account for the 
unexpectedly high winter ISCCP doud amounts. 
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Monthly Mean Total Cloud Amount ISCCP-C2 and WAR 

Ocmt Atm 63 -90 
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Fig. 1: Annual cycle of total cloudiness in the Arctic (north of 62.5° N latitude) from the 

ISCCP-C2 and WAR data sets. ISCCP data are for 1983-1990; WAR data cover the 
period 1951-1981. 


1.2 Radiative Fluxes Estimated from the ISCCP Cloud Data 

Information on radiative fluxes in the Arctic is sparse and originates mainly from a small number 
of ice islands coastal stations. More detailed information with significant spatial coverage cannot 
be expected within the next 10-20 years. In order to quantify the Arctic radiation balance for 
large-scale process studies such as sea ice modeling experiments as well as for comparison with 
GCM modeling studies, we must rely on satellite remote sensing techniques which provide the 
only source of data with sufficient spatial and temporal coverage. A data set of monthly Arctic 
Ocean radiative fluxes at the surface and top of the atmosphere has been compiled for the years 
1983-1986 (Schweiger and Key, 1992). This "climatology” will provide an alternative to prior 
estimates and incorporates the most recent satellite-derived data sets. 

Radiative fluxes are computed using the ISCCP-C2 data set as the main source of information for 
radiatively important atmospheric and surface properties. The ISCCP-C2 data set is a 
compilation of monthly statistics from the 3-hourly ISCCP Cl data set and contains information 
on satellite-derived cloud fraction, atmospheric profiles of temperature, water vapor and ozone, 
surface temperature and reflectivity, cloud top pressure and temperature, and cloud optical depth. 
The compilation of this data set is an ongoing effort, but only data covering the period July 1983 
to December 1986 are used in this study. This data set provides global coverage from a variety 
of sensors on a 250 by 250 km analysis grid. For the polar regions the data set consists entirely 
of data from the AVHRR and TOVS instruments on board the NOAA satellites. The compilation 
of this data set involves a large number of processing steps; for more detail the reader is referred 
to Rossow et al. (1988). 

The radiative transfer model used to calculate radiative fluxes in this study was modified from 
Tsay et al. (1988) by replacing the discrete ordinate solution of the radiative transfer equation 
with a delta-Eddington two-stream approximation for short wave calculations and a hemispheric 
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Shortwave Flux at the Surface 




Fig. 2: Annual cycle of shortwave fluxes 
(Wm -2 ) for ISCCP baseline case. 


Fig. 3: Annual cycle of longwave fluxes 

(Wm* 2 ) for the ISCCP baseline case. 


The spatial variability was also examined (Fig. 4 and 5). In March the central Arctic regions 
receive increasing southward up to 150 Wm* 2 in the Bering Sea area. A modulation of the solar 
zenith angle effect by cloud amount is apparent in the Greenland/Norwegian Seas area, where 
cloudiness is usually greater than in other regions of the Arctic. In April this pattern persists but 
downwelling fluxes have increased by about 100 Wm 2 for the areas north of 70 with a 
latitudinal gradient of 100 Wm' 2 similar to the one found in March. By May, when incoming 
solar radiation at the top of the atmosphere varies little with latitude, cloud amount becomes the 
dominant factor determining downwelling shortwave radiation at the surface. Since the ISCCP 
data set reports fewer clouds in the central Arctic than elsewhere, this region receives more solar 
radiation (300 Wm' 2 ) than other areas of the Arctic. The minimum again occurs in the 
Greenland/Norwegian Seas area where cloud amounts are persistently high. In June and July, 
downwelling shortwave patterns vary little over most of the Arctic region with downwelling 
fluxes in the order of 300-325 Wm' 2 and 275 Wm 2 respectively. A steep gradient can again be 
observed in the Greenland/Norwegian as well as the Barents Seas, where summer cloud amounts 
of 80% reduce incoming shortwave radiation substantially. In these areas summer ice 
concentrations, and therefore aerially averaged albedo, are lower. Since clouds reduce 
downwelling shortwave radiation more effectively over lower albedo surfaces, the combination 
of greater cloud amounts with lower surface albedo leads to downwelling shortwave fluxes as 
low as 150 and 125 Wm' 2 in June and July over the Norwegian sea. In August, the solar zenith 
angle-related concentric pattern redevelops and cloudiness has a secondary effect on the spatial 
distribution of shortwave fluxes. 


Downwelling longwave fluxes show very little variation over most of the Arctic and are on the 
order of 175 Wm" 2 from November through April. A strong increase up to 300 Wm can be 
found over the Norwegian, Greenland and Barents Seas areas. This gradient, following roughly 
the Greenwich meridian is also fairly constant for these month and is again caused by the greater 
cloudiness over these areas. In May the peripheral seas, because of rising temperatures, begin to 
receive greater amounts of longwave radiation (200 Wm' 2 ), while the central Arctic remains at 
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the 175 Wrrr 2 level. In June a sudden jump to 250 Wm -2 with little variation over the central 
areas occurs, and the maximum is reached in July with 275 Wm* 2 over most of the Arctic Ocean 
and the peripheral seas. In August this pattern persists at approximately the same level, before 
values drop again, most drastically in the central Arctic. 

Despite sig nifi cant uncertainties in accuracy, the ISCCP data set for the first time provides an 
opportunity to calculate radiative surface fluxes and cloud forcings from satellite at a spatial and 
temporal resolution suitable for sea ice modeling experiments. The accuracy and potential errors 
of the calculated fluxes are under investigation through comparisons with surface measurements 
and sensitivity studies. 


Downwelling at ® le ® ur * ace 



Fig 4: Spatial variability of downwelling shortwave fluxes (Win' 2 ) for June 1984 - 1990 as 
derived from ISCCP 


OownweHiufl *• 



Downwelling Lon^ave^FJjwgs at the Surface 



Fig 5: Spatial variability of downwelling longwave fluxes (Wm* 2 ) for January and June 1984 - 
1990 as derived from ISCCP 
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1.3. Radiative Fluxes and Cloud Forcing 

Information on radiative fluxes in the Arctic is sparse and originates mainly from a small number 
of ice islands and coastal stations. The monthly cloud data product of the International Cloud 
Climatology Project (ISCCP) has been used to estimate Arctic radiative fluxes and cloud forcing at 
the surface and top of the atmosphere. An example of the results is shown in Figure 6. which gives 
the shortwave, longwave, and net cloud forcing at the surface. Throughout much of the year the 
presence of cloud covers tends to warm the surface. During summer, however, clouds significantly 
decrease the downwelling shortwave radiation and cool the surface. Despite significant 
uncertainties in accuracy, the ISCCP data set for the first time provides an opportunity to calculate 
radiative surface fluxes and cloud forcings from satellite at a spatial and temporal resolution 
suitable for sea ice modeling experiments. The accuracy and potential errors of the calculated 
fluxes are under investigation through sensitivity studies and by comparisons with surface 
measurements. 


Cloud Forcing at the Surface 



Fig. 6: Shortwave, longwave, and net cloud forcing at the surface in the Arctic. Data are 

computed from the ISCCP monthly cloud product for the years 1983- 1990. 
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2. SNOW AND ICE ALBEDO 


Konrad Steffen, Axel Schweiger and Jeffries Key 

2.1 Calibration of AVHRR Channel 1 and 2 

2.1.1 Objectives 

To discuss the Ration (AVHRR) channel 1 and 

"cSr^ocedureu^d ^ound based a~ph«dc 

m estimate *e Sifidence level for the parameterization of surface energy 
flux components as derived from satellite data . 

2.1.2 AVHRR calibration 
NOAA-ll AVHRR LAC data (level 

format includes navigation and mt rellK)ved by linearly nndging the 
Small registration ««••£•-> -*• “ errors can men be reduced to ±1/2 

::r r ■* *. da* - 

Lear ftmcdon: a scan angle of 30 degrees resnlts in a pixel resolution of .5 km. 

2.1.3 Surface and planetary narrow-band albedo 

_ ... . - the avhRR imagery is channel dependent. For channel 1 and 2 pre-launch values 

<DN> to percent albedo (Brown. 1988 k For each pixel tn 
SS were computed and cosine correcnoos were applied. 

K. planetary nanow-baM alfcdo for AVHRR ^ on"L 
were compared against ground based 1 155 m a s 1 ). An integration sphere 

Greenland ice t « ^ ^c —ems. directeTwanls the 

was mounted in front of the optical eaa v , reflectance. The snow 

aemth and pe^dicnla, n> The prevailing air 

was less than one day old. wnth gram broadband global radiadon, diffuse and reflected 

temperature was well below freezing, in addition, broadband glow ™ ^ radiation 

shortwave radiadon - *L? hom me 

balance expenmen t at me Legrated reflectance for the same spectral region was 

STS. mm broadbLd albedo nre^nremems from pyranomerer and 

spectroradiometer can be achieved within 0.7 % accuracy. 

The atmospheric absorption 

LOWTRAN 7. This work was performed by Mr. Marcel Haiuger. cm 
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runs were performed: (a) using the in sirn radiosonde ~ 

A« speerrai 

band-width (^-1^2 )* ^ following. 

| [S(X)<t>i(X)x(A) a g (A)]dA/J [S(X)4>i(A)]dX 
where S(A) is the solar 

function of the AVHRR channel (i=l,2), a g (A) is the nanw o q ^ of te 

(X) is the calculated and the satellite. The integration was 

atmosphere and the surface, and ( ) .. . nrofile data and (b) standard sub-arctic 

done for Mb . AVHRR cbarn£ and calorlated 

atmosphere (Tab- „ dL0 WTRAN) narrow-band planetary albedo for channel 1 is 0.5% using 
(spectroradiometer and L °^R) ^ 7 model Us ing the standard sub-arctic atmosphere in 
radiosonde profile data in the LOWTR AVHRR channel 2 the calculated 

the LOWTRAN model, the difference increased to 2.9%. For AVMKK cnauuc 

^ melred narrow-band albedo values differ by approxrmamly 9%. 

Surface narrow-band albedo a g can also be derived from planetary narrow-band albedo a p values 
as following: 

ag = Up / [x (sun-surface) ^ "^(surface-satellite) 

The transmission values were calculated using 

intrinsic atmospbenc reflectance due to e derivod AV HRR planetary 

arctic atmosphere. tZ of magmtude for me crforlated 

^Xo“w— o“n^rprising as txnh —os trie dm LOWTRAN 7 

model. 

Tawe i: 

ground-based spectrometer vafues a g and radiative transfer modebng. 

I AVHRR a g I calculated a p a P I H 

Channel AVHRR with radiosonde atm. with smndmd ann- . S P‘” r ^ g 

l l380dS80nmr~ 68.8% 81.0% 

2 (725 - 1 100 nm) ■ 1% ' 65.2% 1 J 


“P 

AVHRR 

68 . 8 % 

56.1% 


calculated ap 

with standard atm. 
65.5% 
64.5% 


a g 

setrometer 

91.6% 

81.0% 


™' e 2: 

and X2(surface-satellite). 

avhrr [ T % calculated a g calculated a g "'l T 2 

Channel Spectro- with radiosonde withst ^^ atm 

meter AVHRR AVURR - 


AVHRR I 

a g 

Channel 

Spectro- 

meter 

l 1 (580-680 nm) 

91.6% 


81.0% | 


0.825 


0.878 


0.895 


0.922 
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2. 1.4 Conclusions 

Standard sub-arctic aerosol values were used for the computation of aerosol scattering as 
atmospheric profile measurements were not available. This might explain in part the discrepancies 
between measured and modeled surface albedo values. Aerosol content can vary with time and 
space significantly which can reduce the reflectance by several percents. Assuming that the 
radiative transfer model is correct, we can speculate that a major drift in AVHRR channel 2 has 
occurred since the launch of NOAA 1 1 (November, 8, 1988). For the visible channel of NOAA 9 
degrading rates of up to 6 % per year were reported (Staylor, 1990). This case study should be 
judged as a preli min ary analysis, as more cases have to be compared for a statistical analysis. In 
addition, LOWTRAN 7 derived spectral transmission values have to be validated to separate the 
modeling error from sensor drift. 

Narrow-band surface albedo values can be estimated from AVHRR channel 1 and 2 measurements 
with an accuracy of 2%, 14% respectively using LOWTRAN 7 radiative transfer modeling and 
radiosonde profile information. More case studies will be analyzed to discuss this technique in 
more depths and to derive statistically significant values. 


2.2. Spectral Reflectance 
2.2.1 Hemispheric Reflectance 

Hemispheric spectral albedo was measured for five different snow grain size ranges (Table 3, Fig. 
7). The largest spectral reflectance in the visible spectrum was found for very new snow. In the 
near infrared region (900 - 1500 nm) the reflectance dependents on the soot content, a parameter 
not measured in situ during the field experiment. Based on the spectral characteristics of different 
snow surfaces, satellite and aircraft measurements with high spectral resolution have the potential 
for classifying the different snow fade’s. 


Table 3: Snow properties for different snow types given in Figure 7 


Snow Type 

Age 

Grain Size (mm) 

Solar Zenith Angle 

Dry Snow 

2 hours 

0.5 

50 

Dry Snow 

1 day 

0. 1-0.2 

50 

Dry Snow 

3 days 

0. 1-0.3 

50 

Wet Snow 

- 

1.0 

50 

Wet Snow 

- 

2.0-5.0 

68 


2.2.2 Anisotropic Reflectance 

The anisotropic reflectance factor (ARF) was calculated based on the field measurements for dry 
snow with grain sizes varying between 0.1- 0.3 mm (diameter). Figure 8 shows the ASF for a 
reflectance an gle of 85° and a solar zenith angle of 71°. The results are in good agreement with 
earlier measurements taken in the spectral range 400 - 500 nm (Steffen, 1987). The measurements 
show that dry snow can be treated as an Lambertian reflector with minor spectral dependency. The 
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results for wet snow have not fully been analyzed yet, but they clearly show a strong specular 
reflectance dependency on the azimuth angle. The wet snow ARF values will be needed for the 
conversion of satellite-based directional-spectral reflectance to albedo. 



Fig. 7: Hemispheric spectral albedo for very new snow, new snow, onset of melt, and melting 

snow surface in the range 300 - 2500 nm. Snow properties and solar zenith angles for 
each measurement is given in Table 3. 


Fig. 8: 
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Greenland; 23 May 1991 
Solar Zenith: 71 
Reflectance : 85 
Snow cond.: dry 
Snow grain: 0.1 -0.3 mm 
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Anisotropic reflectance factor for seven different azim uth angles, a reflectance angle of 
85°, and a solar zenith able of 71°. The snow was dry with grain radii varying between 
0.1 - 0.3 mm. 
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2.2.3 Directional hemispheric reflectance 

Bidirectional Distribution Reflectance Function (BDRF) was modeled replacing the two-stream 
approach as discussed by Warren and Wiscome (1981) with a discrete-ordinate solution (48 
streams). The model assumes Mie-scattering of optically equivalent ice spheres. Small amounts of 
carbon soot are included in the medium and are needed to achieve a match between measured and 
modeled reflectance's in the visible region. Presently modeling results assume a single homogenous 
layer of snow but multi-layer calculations are also possible. 

The model results were compared with field measurements as shown in Figure 9. In general there 
is a good agreement if a soot amount of 0.30 ppmw were used. This soot amount seems rather 
high for the Greenland region, however, in situ measurements are missing and can not confirm this 
assumption. If the model run was performed with 0.00 ppmw soot amount, the spectral reflectance 
in the visible range (0.3 - 0.8 mm) was approximately 0.1 higher than the in situ measurement. 


1.0 1' ! 1 I 1 ' 'T 1 ' I ' ' 



Wavelength [jimj 

Fig. 9: Directional hemispheric albedo of snow from with discrete-ordinate model versus ground 

measurement for a dry snow surface with a mean grain size of 0. 1 mm. For the model 
computation a soot amount of 0.3 ppmw was assume d 


23 Comparison of Two Satellite-Derived Albedo Data Sets 

Variations in the surface albedo of the Arctic pack ice cover is an important factor in the long- 
term mass balance and stability of the pack ice. Continued monitoring of albedo is consequently 
an important component of the "early detection" strategy for identifying global and regional 
climate changes. Nevertheless, until recently, studies of the Arctic Ocean albedo have been of 
limited spatial and temporal scope, based largely on observations at drifting stations, fast ice, and 
during aircraft missions. Other investigators have used these measurements and satellite passive 
microwave data to estimate regional summer albedo (Robinson et al., 1992). 
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With the aim of developing a long-term record Robinson et al. (1992) recently developed a ten 
year data set of parameterized albedo for the entire Arctic Ocean for May through mid-August. 
Maps of surface brightness charted through manual analysis of visible-band DMSP imagery were 
converted to parameterized surface albedo. Although lacking radiometric control, this data base 
(hereafter referred to as the R092 set) has allowed for definition of the climatological 
characteristics of albedo. More recently, we have compiled a separate albedo data set, based on 
the fluxes derived from the ISCCP monthly data set as described in the previous section. We 
have compared these two data sets for the months of May, June, and July. While providing an 
evaluation of uncertainties in our current knowledge of Arctic Ocean surface albedo, such a 
comparison also makes clear problems in the automated extraction of these data, and clarifies 
aspects of the existing retrieval algorithm in need of further attention. 

Results may be summarized as follows. During May, when snow melt is largely confined to 
coastal regions (Robinson et al., 1992) both data sets show little spatial variability over the 
central Arctic. The R092 albedo over the central Arctic tend to be between 0.75 and 0.80 (the 
latter figure representing the maximum possible parameterized value). By contrast, the ISCCP 
values over the same regions are lower by about 0.05 to 0.10. Both data sets show a strong 
gradient of decreasing albedo near coastal regions, especially toward the Norwegian and East 
Siberian Seas, consistent with the effects of decreasing ice concentrations and coastal snow melt 
(Robinson et al., 1992). Both data sets also show surface albedo to be lower near the New 
Siberian Islands. The reduction in albedo near the western portion of the Canadian Arctic 
Archipelago shown by ISCCP, however, is not depicted in the R092 analysis. 

By June, melt is underway over large areas of the sea ice cover (Robinson et al., 1992) which 
results in a decrease in surface albedo (Fig. 10). Both data sets show a concentric pattern of 
surface albedo, with values decreasing from the pole southwards. Serreze et al. (1991) show that 
this concentric pattern can be related to the climatological distribution of melting degree days. 
ISCCP albedo, however, again tend to be lower than those reported by R092. For example, in 
the central Arctic regions ISCCP tends to be lower by about 0 . 1 . 

Albedo continue to decline during July. The concentric pattern observed in both data sets for 
June is still present. During July the ISCCP and R092 albedo match very well, except near the 
pole where ISCCP albedo are now higher than the R092 values by up to 0.05. 

In summary, generally good agreement in the spatial patterns of surface albedo is found between 
R092 and ISCCP data sets. While the R092 data set cannot be considered "ground truth", the 
similarity between albedo patterns in these two independently-derived products gives a positive 
outlook for operational retrieval of albedo from satellite data. However, there are systematic 
differences in the absolute albedo values, with ISCCP albedo being lower during May-June and 
higher in July. Reasons for these differences presumably reflect the present inadequacies in the 
ISCCP cloud-clearing algorithm and the simplifications and subjectivity in the R092 procedures. 
Because of these uncertainties it is not currently possible to determine which is more correct. 
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3. LEADS AND MODELING 

Jim Maslanik and Jeffries Key 

3.1 Effect of Measurement Scale on Turbulent Fluxes from Leads 

Fractures in the sea ice pack ("leads"), either open or refrozen, are an important component of 
local scale heat exchange in the Arctic, providing a significant source of heat and moisture to the 
atmosphere that is as much as two orders of magnitude greater than turbulent energy exchange 
from thick sea ice. Lead coverage also affects the ice-pack albedo and in turn, the rate of ice melt 
through the absorption of heat within the open-water leads (Maykut and Perovich, 1987). 
Modeling studies of the influence of sea ice leads on climate point out the importance of 
accurately representing lead fraction, which typically defines the proportion of open water within 
the interior pack ice (Ledley, 1988). Relatively large changes in lead coverage can be related to 
strong storms in the Arctic (Maslanik and Barry, 1989) and can contribute as a positive feedback 
to maintaining the strength of these storms. In addition to total lead-covered area, the widths of 
individual leads affects the rate of heat transfer from the ocean to the atmosphere (Andreas and 
Murphy, 1986). Unfortunately, relatively little information is available on the numbers, 
dimensions, and life cycle of leads in the Arctic ice pack. Remote sensing - in particular the use 
of medium-resolution imagery from polar orbiters such as the NOAA Advanced Very High 
Resolution Radiometer (AVHRR) -offers the potential to develop a climatology of lead statistics 
that is well-suited for climate studies and comparison to other remotely-sensed data such as 
microwave imagery. However, the validity of the retrieved geometrical characteristics of leads - 
in particular lead widths and area fraction - from medium and low-resolution data is uncertain, so 
that any estimates of their climatic impact are of questionable utility. 

Key et al. (in press) have examined the effects of satellite sensor resolution on lead width 
distributions, mean lead widths, and total area coverage of leads within an image. It was shown 
that small leads tend to disappear in the coarser resolution data and large leads "grow" primarily 
by becoming wider. The total lead fraction decreases as field-of-view increases, roughly 
following a log-linear relationship. Both the rates and direction of change are, however, closely 
tied to the threshold used in creating the binary image. Lead orientation distributions also change 
with increasing pixel size indicating that, in the imagery examined, small leads exhibit different 
orientations than larger leads. 

How do these relationships affect area-averaged estimates of turbulent heat flux? Changes in 
both the mean lead width and lead-covered area are considered in the calculation of sensible and 
latent heat flux as a function of fetch (treated here as the lead width), surface temperature, air 
temperature, and wind speed using the procedure outlined by Andreas and Murphy (1986). In 
this approach, a bulk Richardson number defines atmospheric stability that controls convective 
turbulence based on temperature and wind speed. Convective turbulence combines with the 
mechanical mixing introduced by the step effect of an air mass in equilibrium with thick sea-ice 
conditions. This air mass is traveling over the physically rough edge of a lead and the 
considerably warmer open water or thin ice in the lead. The addition of mechanical turbulence 
introduced by the ice-lead boundary tends to result in a higher rate of heat transfer from smaller 
leads compared to larger leads. Thus, for a given aerial coverage of leads in an image, a greater 
number of smaller leads will result in more heat loss to the atmosphere than from a lesser number 
of larger leads, even though the total amount of open water in the image remains the same. 
Under the conditions examined by Andreas and Murphy (1986), this decrease in flux as lead 
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width increases becomes negligible for lead widths greater than about 200 m. 


To illustrate the effects of changes in lead statistics using different image fields-of-view we 
calculated sensible and latent heat flux using the above approach with a variety of Landsat 
images successively degraded from a pixel size of 80 m to 160, 320, and 640 m pixel size 
images. An open-water temperature of -1.8° C, wind speed of 5 m s -l. air temperatures of -28 9° 
C at a reference height of 2 m, ocean salinity of 34 ppt, air pressure of 1000 mb and a neutral- 
stability drag coefficient of 1.49 x 10'3 are used to represent typical mid-winter oZJv) 
conditions over the Arctic sea ice pack. Although leads are often cohered by thin ice rather tiUn 

2*? ™“ er J us have * lower surface temperature than open water, the assumption that the 
eads are not refrozen and have a surface temperature of -1.8° C is a useful baseline for our 

(sensible P lus latent > heat flux from leads is calculated using the mean 

3t 6aC n ° f ViCW md then wei 8 hted b y ^ aerial coverage of leads for the images 
to yield an aenally-averaged heat flux. Turbulent fluxes from open-water leads under these 
conditions are around 300 Wnr* compared to a flux of nearly 0 froiTsurrounding ice taken To te 
three meters thick. Thus, lead fraction and lead width dominate the transfer o § f turbulent he^ 
fltrough the ice pack during winter. Table 4 shows these aerial averages for six MSS images 
Since the effect of increasing the fields-of-view in these examples is to decrease the apparent tead 
fraction, aenafly-averaged fluxes decrease as field-of-view increases. However, as nofed earlief 
the choice of tlvesholds can affects both the magnitude and direction of change in lead statistics 

FOV ' ! . f we that the lead w »dths and lead fractions measured using the 80 

m FOV imagery are closest to reality, then the errors introduced by using lead widths ^rd lead 

:r n “, f 3 m F ° V “ SUbSt3ntial - avera ^ 45% over ^he six images SiS 
tte^anse m turbulent heat transfer with changing lead width is greatest for smaller leads those 

37? S “ aUeSt mean lead widths at the 80 m FOV (such as Images C and F) are most 
affected In die images studied here, where the mean lead width is fairly large the effect of 
errors in lead fraction is about five times that of the effect of uncertainty in iLd wfdth 


Table 4. 


Aerially averaged turbulent (sensible and latent) flux (in Wm*^ fnr rvni^i i an a- • 

view forsiY MSS , ) x {m wm ) tor typical January conditions as a function of field of 

MSS images. The percent change in the flux between FOVs of 80 m and 640 m is also showu. 


JFOV (m) 

A 

B 

C 

D 

E 

F 

80 

23.5 

19.1 

17.4 

50.7 

10.6 

15 0 

160 

22.2 

17.7 

14.4 

47.0 

9.7 

12 1 

320 

19.9 

16.1 

11.7 

41.6 

8.2 

8.9 

640 

16.7 

14.0 

7.9 

33.7 

5.6 

5.3 

% change 

29 

27 

55 

34 

47 

65 


3.2 Modeling 

^ - V 7j Uable Way ° f assessin 8 ** relati ve importance of remotely-sensed parameters and of 
assigning accuracy requirements to remote-sensing algorithms is through numerical sensitivity 

of ic^e°ad t tl!m^ n ° SPher 1 SyStCm ' ^ 6Xample ° f studies > we focus 00 defining the rok 
ead geometry in simulations of the sea ice cover. Low-resolution imagery such as AVHRR 

nor nifiTf mcr0wave data provide formation on total open-water area within regions but are 
not useful for mapping the widths of individual leads. Leads can be mapped in hSh resolution 

^of^T? “ h AR ’ bUt , 0nly With cons “ teah,e effort. However, research has shwvn that the 
rate of turbulent heat transfer varies significantly with lead width for leads less than about 100 m 
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The question arises, then, as to whether accurate simulations of sea ice require lead width data. To 
test this and other such sensitivities, we use a sea ice model that simulates the Arctic ice pack, and 
includes a coupling to the atmospheric boundary layer. This coupling is important in that it allows 
changes in the ice pack to affect the overlying atmosphere, whereas previous sensitivity studies 
typically do not allow this feedback. The ice model is run using daily wind fields, climatological 
ocean currents, snowfall, and cloud properties, and modeled longwave and shortwave radiation. 
Typically, our model runs for this project use daily forcings for 1989-1991. The effect of lead 
width is addressed by using a larger bulk heat transfer coefficient for narrow leads. 


Results of ice-pack simulations addressing the question of the importance of lead width are 
summarized in Figure 11. The significance of treating lead width in the model depends on the 
desired output. Ice growth remains steady as lead width changes in this approach, due to a 
negative feedback between turbulent heat flux and air temperature. Results show that, when heat 
from leads is removed from the system (e.g., not allowed to warm the near-surface air), ice volume 
increases by 17% for the period showa However, when heat rising from the leads is retained in 
the boundary layer, the effect of lead width (or the effect of increasing the rate of heat transfer) is 
3%. This difference is due to the fact that, although the heat transfer is more efficient from small 
leads, surface air temperatures are higher, so the rate of heat loss and the resulting ice growth rate, 
are relatively unaffected. Thus, in terms of ice production in this type of model, the importance of 
lead width lies somewhere between these two extremes where all heat is removed from the system 
versus the cases where heat is retained in the atmosphere above leads and sea ice, with no lateral 
mixing of air ova the two surfaces. In terms of atmospheric effects and climate modeling, taking 
into account the more efficient heat transfer from narrow leads produced mean increases in surface 
air temperature of about 0. 1 C in the central Arctic and approximately 1.2° C in the seasonal sea 
ice zones in the extreme case tested here of no air-mass mixing or advectiom This difference is due 
to the proportion of open water in the ice pack, which averaged less than 2% in the central Arctic 
in these simulations. In particular, these results suggest that in terms of turbulent heat transfer 
(ignoring effects of lead width on solar heat absorption and lateral melt rates), detailed knowledge 
of lead width distributions may be of secondary importance. More such experiments are planned 
to test the sensitivity of the Arctic ice-atmosphere system to other remotely-sensed parameters, as 
well as to a more realistic situation where lateral mixing of air over leads and ice is permitted. 


Fig. 11: A comparison of the effect of 
including lead width (taken as equal to 
fetch) in the calculation of turbulent heat 
flux in simulations of the Arctic ice cover. 
Lead width parameterized as a function of 
ice concentration, with heat transferred 
more efficiency from narrow leads than 
wide leads. Ice volume represents the 
sum of all ice in the Arctic and peripheral 
seas. Simulated ice volumes are shown 
for the period from 1 Jan. 1989 - 3 Nov. 
1990 to represent the seasonal cycle. In 
the "No ABL" cases, the coupling between 
the sea-ice surface and the overlying 
atmospheric boundary layer is turned off, 
so that surface air temperatures are 
independent of surface heat flux. In the 
M ABL M cases, the coupling allows heat 
from the surface to warm the overlying air 
to separate equilibrium temperatures over 
ice and open water. 


Effect of Fetch on Ice Volume 
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